Dynamics of wrinkles on a vesicle in external flow 
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Recent experiments by Kantsler et. al. (2007) have shown that the relaxational dynamics of a 
vesicle in external elongation flow is accompanied by the formation of wrinkles on a membrane. Mo- 
tivated by these experiments we present a theory describing the dynamics of a wrinkled membrane. 
Formation of wrinkles is related to the dynamical instability induced by negative surface tension 
of the membrane. For quasi-spherical vesicles we perform analytical study of the wrinkle structure 
dynamics. We derive the expression for the instability threshold and identify three stages of the 
dynamics. The scaling laws for the temporal evolution of wrinkling wavelength and surface tension 
are established and confirmed numerically. 
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Wrinkling of different thin sheets is a weU-known effect 
which can be frequently observed in everyday hfe. Usu- 
ally wrinkle patterns appear due to the external tensions 
applied to some material, or as a result of compression 
of inextensible films. Main properties of steady and/or 
equilibrium wrinkling structures are now well understood 
[if, Q ■ However, much less is known about the dynam- 
ics of wrinkle structures. In certain biologically moti- 
vated experiments on membranes and vesicles 0, , the 
wrinkles are formed because of the instability induced by 
negative membrane tension, which is closely related to 
a buckling instability. Wrinkles formed this way exhibit 
non-trivial growth and relaxational dynamics. Theoret- 
ical description of these essentially non-equilibrium pro- 
cesses is a challenging problem which we attempt to ap- 
proach in this paper. Although in this letter we focus on 
the analysis of the recent experiments 0], some predic- 
tions of our theory are universal and may be successfully 
applied to other systems where the formation of wrinkles 
is caused by negative tension. 

Vesicles exhibit a variety of different regimes of motion 
in external fluid flows. These regimes were extensively 
studied during the last decade bot h by the experimental- 
ists [1,11, 0,i and theoreticians [1 [M (H!, E 111 [13, [ll ■ 
Depending on the external parameters a vesicle in exter- 
nal shear flow can experience three types of behavior: 
tank-treading, tumbling and trembling (also referred as 
vacillating-breathing 13] or swinging [l2|) [sj. Although 
these motions correspond to c^uite non-trivial dynamics, 
the shape of a vesicle always remains smooth and can be 
effectively approximated by an ellipsoid. Recently, ex- 
periments performed by Kantsler et. al. [^ revealed a 
qualitatively new effect observed in non-stationary elon- 
gation flows. It was shown that in strong flows the re- 
laxational dynamics of a vesicle is accompanied by the 
excitement of high order membrane deformation modes 
called wrinkles. Such dynamics can not be described in 
a framework of low-dimensional models which were used 
for the analysis of vesicle tank-treading, tumbling, and 
trembling. In this letter we extend these models to in- 
clude the interaction between the vesicle shape and the 
wrinkle structure. 



Before proceeding further we would also like to men- 
tion the experimental flG*! and theoretical [l3| investiga- 
tions of wrinkle formation on microcapsules in external 
shear flows. Although this effect is similar to the one 
discussed here the underlying physics and main proper- 
ties of wrinkles are essentially different. For example, the 
wrinkles which are observed on vesicles are not stationary 
and are excited only for a limited amount of time. 

This letter is organized as follows: First we discuss 
the main features of a vesicle and show that the nega- 
tive tension leads to instabilities of the flat membrane. 
Second, we present a model of a quasi-spherical vesicle 
in external flow. We derive the threshold of the insta- 
bility and analyze the dynamics of wrinkle formation in 
strong flows. We show that one can distinguish three dif- 
ferent stages of the dynamics. The first one lasts for a 
vanishingly small amount of time and is characterized by 
the exponential growth of short wavelength excitations. 
During the second stage the amplitude of wrinkles and 
their characteristic wavelengths grow algebraically. We 
derive the corresponding exponents and compare them 
to the experimental results. During the third stage the 
amplitude of wrinkles gradually reaches maximum and 
afterwards decays to zero. We find the scaling estima- 
tions for the wrinkle wavelength and the surface tension. 
In the end we present the results of our numerical simu- 
lations which confirm the analytical theory. We conclude 
by discussing the future challenges. 

Vesicles are closed lipid bilayers which are incompress- 
ible and impermeable to the surrounding fiuid. These two 
properties result in the conservation of the vesicle volume 
and the membrane area. Thus any vesicle is character- 
ized by its excess area A which is the measure of vesicle's 
"nonsphericity" : A = {An + A)R^ where A is the mem- 
brane area and R is the effective vesicle radius, defined 
by the vesicle volume: V = AttR'^ /3. Free energy of the 
closed membrane is defined by the Helfrich functional and 
consists of the contributions from the bending energy and 
the surface tension a [Tsj : 

F^JdA[^H' + a). (1) 
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Here H is the local mean curvature and k is the bending 
rigidity of a membrane. Note that for closed membrane 
geometry, the surface tension cr is a quantity adjusting 
to other membrane parameters (similar to the pressure 
in an incompressible fluid) to ensure a given value of the 
membrane area A. It is useful to analyze the stability of 
a flat membrane with a given value of a before proceed- 
ing to the case of a closed membrane with a fixed area. 
Small perturbations of a flat membrane can be param- 
eterized by a height function z = h{x,y) which can be 
expanded in Fourier harmonics: h(r) = ^ /ik exp(ikr). 
The quadratic part of the Helfrich energy has the follow- 
ing form: 



(2) 



For positive cr > 0, this function is minimized by m = 0, so 
the flat membrane state is stable. However, for negative 
tension a < 0, when the membrane is being shrinked, 
the modes with k < \/\<j\/k become unstable. As we 
will show, this particular kind of instability is responsi- 
ble for the formation of wrinkles in the experiment 
More precisely, we will show that the surface tension be- 
comes negative when the direction of the external flow is 
reversed instantly. In this case small thermally induced 
membrane deformations on top of a smooth vesicle shape 
are amplified and lead to the formation of wrinkles. 

In order to study this effect quantitatively we ana- 
lyze the model of quasi-spherical (A ^ 1) vesicles which 
proved itself to be very successful in analytical investiga- 
tions of vesicle dynamics in external flows [3, [3 ■ 
The vesicle shape is parameterized by the small displace- 
ment function u{6,4'): r = R{1 + u) which can be ex- 
panded in the series of spherical harmonics: 
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UlUt)yim{0,(f') (3) 



The external velocity is assumed to have a linear elon- 
gational profile with the time-dependent strain: dxVy = 
dyV^ = 11^5/(1676^) • S{t)VA/T, where t = r]R^/K 
is the characteristic time-scale associated with the mem- 
brane bending forces. The numerical factor in this def- 
inition was included to simplify the expressions below. 
Throughout the paper we will discuss the experimentally 
interesting situation when S{t) is the Heaviside-like func- 
tion S{t) = -S'sign(i). 

The dynamical equations describing the dynamics of a 
quasi-spherical vesicle in external flow were first derived 
in in the leading order in the small parameter y/A <^ 
1: 

TUirn = S{t)fi„i - {Aia + Ti)uirn + Om(i), (4) 

Here fim — Si^2iSm,2 + ^^,-2) and a{t) is the dimcn- 
sionless angularly averaged part of the surface tension 
which is a Lagrangian multiplier associated with the ex- 
cess area conservation constraint ^ |wimP = 1- The 



numerical coefficients of ^ are given by: Ti ~ {I — 
+ 1)2(/ + 2)/{2l + 1)(2P + 21 - 1), and Ai = 
l{l + + / - 2)/(2/ + l)(2/2 + 21 - 1). The statis- 
tical properties of thermal Langevin forces (t) can be 
found in [9]. In this letter we assume that the tempera- 
ture T is small, so that Ci,n{t) in (HI are negligible dur- 
ing the evolution of the wrinkle structure. They only 
affect the dynamics through the non-zero initial condi- 
tions u;.m(t = 0) which ensure that the instability can 
indeed take the system out of the unstable state. In 
stationary planar elongational flows most of the excess 
area is stored in the 3^2, ±2 harmonics. We therefore use 
the following parametrization : U22 — U2 -2 — ^ and 
uim — \/l — |C/p wim for l,m ^ 2, ±2. The fraction 
of total excess area stored in the second-order spherical 
harmonics 3^2±2 is given by \U\'^. The argument of U is 
related to the vesicle orientation angle $ in the xy-plane: 
U = |?7| exp(— 2i$). The normalization condition can 
be rewritten as ^ l^imP = 1- The dynamical equations 
acquire the following form: 



TU = S{t)-{A2a + T2)U 



(5) 



Using the excess area conservation law one can also find 
the expression for the surface tension: 



S{t)RcU ^T2\U\^ -fjl - \U\^) 
A2\U\^ + A{1-\U\'') 



(6) 



where A = J^M^iml'^ and F = J2^i\^im\'^ ■ Using the 
expression ([6|) one can easily show that there will be an 
instability for large enough values of S. Indeed, for con- 
stant positive S{t) — S a.t t < the vesicle exhibits small 
thermal fluctuations near the stable point U = 1. How- 
ever as S{t) changes sign this state becomes unstable and 
the vesicle starts rotating to the new stable point which 
corresponds to ?7 = — 1. The stability of the membrane 
can be studied by analyzing the expression ([6]). The sur- 
face tension instantly becomes negative after changing 
the velocity field: a{t = -|-0) = - (F2 + S) /A2. Negative 
a can destabilize the high-order harmonics. The explicit 
condition can be found from fj]). All harmonics of order 
up to I become unstable ii Aia + Ti < 0, which yields 
S > Si = A2T1/A1 — T2- For large I ^ 1 one can use the 
expansions Ai - //4,r, - 1^/4, to obtain Si ~ A2P > 1. 
The most unstable mode can be found by maximizing the 
growth increment — F/ — Aia. One obtains Iq = \J 5/3^2 
for strong fiows with S" ^ 1. 

Below the instability threshold, for relatively small val- 
ues of S higher order harmonics are not excited and the 
dynamics of a vesicle can be well described in terms of a 
single V variable. The conservation of excess area implies 
that the dynamics is purely rotational and \U\^ = 1. The 
characteristic time-scale associated with the rotational 
dynamics is estimated as t/S. 

Describing the dynamics of the vesicle above the insta- 
bility threshold is a considerably more difficult problem 
which requires an analysis of a complex nonlinear sys- 
tem Q with large number of degrees of freedom. Fortu- 
nately, it is possible to approach this problem analytically 
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FIG. 1: Temporal dynamics of |?7p, A19 (left grid) and sur- 
face tension a (right grid). Dotted line corresponds to the 
theoretical prediction ([8}. 



for strong flows which correspond to /q ^ 1. Although 
quantitative predictions of such an analysis can be ap- 
plied only to strong flows, the qualitative picture is in 
most aspects the same even for moderate strains S >1. 

For large values of S one can distinguish between sev- 
eral different stages in the wrinkle evolution. At the first 
stage most of the excess area is stored in the second or- 
der angular harmonics and the surface tension can be 
assumed to be constant a — —25/^2- Unstable higher 
order harmonics grow exponentially and the distribution 
of wim becomes centered near the most unstable modes 
with I — Iq ^ VS", so the terms A, f can be estimated 
as A ^ ^/S and f ~ S^^^. This rapid growth saturates 
when the total excess area, stored in higher-order har- 
monics, becomes large enough so that their contribution 
to the surface tension becomes comparable to the con- 
tribution from the external flow. Formal condition can 
be found from (1 - \U\^) ~ S-^^^ < 1. Note that 
at this time most of the excess area is still stored in the 
second-order harmonics. The characteristic duration of 
this stage can be estimated as tS~^^'^. This time is much 
less than the characteristic relaxation time of the variable 
U, which can be estimated as tS^^ from the equation ^ 
as in the situation below the instability threshold. 

Second stage starts at i > tS~^^^ after the expo- 
nential growth has saturated. During the second stage 
higher order harmonics wim provide the main contribu- 
tion to the surface tension ^ which can be approxi- 
mated as (T — —(r/A). Note that it depends only on 
the distribution of excess area over the higher-order har- 
monics wirn- In order to find this distribution for large 
I :s> 1 one can use the formal solution of (|3]): wim{t) — 
CwimiO) exp (—Tit/r + Aip) where C is the normaliza- 
tion constant and P ^ — Jq dt'(j{t')/T. The distribution 
of excess area A; = (l—|J7p) is a narrow func- 

tion of I centered around some which is determined as 
maximum oiAip-Tit/r: U = y/pr/it - ^/S. The char- 
acteristic width of the distribution can be also estimated 
in a usual way: 51 ~ {l^t/r)^'^/'^ . It is much less than 
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FIG. 2: I, versus the strain S. Linear curve corresponds to 
the scaling Z, oc 5'^''^. Inset: Distribution of A; as a function 
of I at the moments tS/r = .7 : 1.3 : 2.2 (right to left). 

for t ^ rS'"'^/^. Narrowness of the distribution allows 
one to find the exact expression for the surface tension: 
(T = —T/A = —11 = —pr/St. Using the definition of p we 
obtain the following closed equation: 

P - p/3i (7) 

One can find its solution using the initial condition for 
the surface tension: a ~ ^/S at t ^ tS^^^"^. This 
yields p = c(i/r)^/'^, where c is a constant of order 
unity. Similarly we obtain a = -(c/3)(t/T)-2/3 ^^^^ 

h = \/c/3{t/T)~^^'^ . Therefore, second stage of the dy- 
namics is characterized by the algebraic decay of the sur- 
face tension and by the narrow spectral distribution A; 
of the wrinkling structure, whose peak smoothly drifts 
towards small /. As the absolute value of Z* goes down 
the external velocity contribution S{t) Re U to the sur- 
face tension in ([6]) becomes important again. Comparing 
different contributions in ([6|) one can estimate the dura- 
tion of the second stage as t/S. In the end of the second 
stage the peak of the A; distribution is centered near 
Z* ^ 5^/"^. This scaling law relates wavelength of wrin- 
kles to the strain A = RS~^^^ and is one of the main 
results of this letter. Note that this wavelength is much 
smaller compared to the wavelength of the initially most 
unstable mode which is given by RS~^^^ . Our scaling is 
in reasonable agreement with the experimental results of 
Kantsler et. al. 

Second stage is followed by third one during which the 
surface tension in ^ is determined both by the exter- 
nal velocity and the wrinkles. The characteristic value 
of the surface tension can be estimated as S"^/^ which 
is much smaller than its initial value (t(0) ~ S. Due 
to a significant decrease of the surface tension, the dy- 
namics of the second-order harmonics during the second 
stage is determined solely by the external velocity term 
S{t) in (O. Therefore one can find the solution of ([5]) 
in the leading order: U{t) = 1 — St/T. The characteris- 
tic amplitude of wrinkles obeys the following simple law: 
^1 - jj/p = ^(2 - St/T)St/T. In order to analyze the 
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evolution of the surface tension a and the wrinkles wave- 
length A ~ i?//, one has to solve the equation rp = a 
keeping all the contributions to the surface tension. With 
the paramctrization p{t) — S~^/'^ f^St/r) the problem is 
reduced to the ordinary differential equation: 



df{x) _ fix) 
dx 



1 



3a; 



2 — a; y xf{x) 



(8) 



with the initial condition x '^/'^ f{x) = c ior x — S ^1"^ 
determined by the second stage. It is possible to solve (|S]) 

2/3 



analytically: /(x) 



logCv^a;(2-a;) 



where 



C is a constant of order unity. Note, that this solution 
adds only logarithmical factors of order log S to the previ- 
ously derived scahngs ^ S^^^ and a ^ S"^/^. Although 
the second stage is formally present in the solution of 
([6|) and corresponds to the limit a; <C 1 we purposely 
separate it from the third one because of its universality: 
during the second stage the surface tension is determined 
solely by the wrinkle structure. Therefore one may ex- 
pect that the scaling laws esc and a oc t'"^!^ will 
also hold for other problems where the initial negative 
surface tension was caused by different means. 

We have tested our results by direct numerical simu- 
lations of equation The total number of harmonics 
in our simulations was Imax ~ 1\fS. The temperature 
which determines the power of the thermal noises in 
dH) was taken to be T = 10~^k. On FIG. [1] one can see 
the results of simulations for S — 10'^. Three temporal 
stages with qualitatively different dynamics are clearly 
distinguishable. The behavior of |J7p is very close to 
quadratic and the dynamics of <j{€) may be well fitted 
by the solution of ([8]). The scaling law oc S^^'^ was 
confimed by simulations at different S (FIG. [5]) 

In conclusion we compile a list of main results pre- 
sented in this paper. Motivated by recent experiments 
we studied the relaxational dynamics of a vesicle in 



an elongational flow. We showed that high order mem- 
brane deformation modes are excited by the negative sur- 
face tension induced by external flow. For quasi-spherical 
vesicles we have found an analytical expression for the 
instability threshold and analyzed the evolution of the 
wrinkle structure. We identified three stages of the dy- 
namics. The first stage corresponds to t < tS~^^^ and is 
characterized by the exponential growth of unstable high 
order harmonics with the characteristic scales of order 
A ~ i?S'~^/^. This rapid growth quickly saturates and 
is followed by the second stage. For rS'"'^/^ ^ t ^ t/S 
the surface tension decays algebraically as a{t) oc t^^/-^ 
and the characteristic wavelength of wrinkles grows as 
A cx i^/'^. Characteristic amplitude of the wrinkles grows 
as y/i. During the third stage which ends at t — 2t/S 
wrinkle amplitude behaves like \/ (2 — St/T)St/T and 
the characteristic wavelength can be estimated as A ~ 

Finally we would like to note, that the algebraic growth 
of the characteristic wrinkle wavelength although with 
different exponent \{t) cx i^/"* was observed in the stud- 
ies of the thin rods dynamics experiencing the buckling 
instability [IBi]- We believe that the algebraic growth 
of wrinkle wavelength is a universal property of the in- 
terface dynamics experiencing the buckling instability. 
Although in this letter we considered a particular exper- 
iment, we believe that the scaling laws for the wrinkle 
wavelength \{t) and for the surface tension a{t) derived 
for the second stage are universal and can be observed in 
other experiments where the total area of the interface 
is conserved and the wrinkle growth is initiated by the 
external forces leading to negative surface tension. 

We are indebted to V. Kantsler and V. Steinberg for 
turning our attention to this problem. We also appreciate 
the fruitful discussions with V. Lebedev. This work was 
supported by the RFBR grant and by Dynasty and RSSF 
foundations. 
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